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Abstract 

We present 2D hydrodynamic simulations of circumstellar 
disks around protostars using a 'Piecewise Parabolic Method' 
(PPM) code. We include a point mass embedded within the 
disk and follow the migration of that point mass through the 
disk. Companions with masses Mc ^ 0.5Mj can open a gap 
in the disk sufficient to halt rapid migration through the disk. 
Lower mass companions open gaps, but migration continues 
because sufficient disk mass remains close to the disk to exert 
large tidal torques. We find that the torques which dominate 
the migration of low mass planets originate within a radial 
region within 1-2 Hill radii of the planet's orbit radius, a dis- 
tance smaller than the thickness of the disk. We conclude 
that a very high resolution 3D treatment will be required to 
adequately describe the planet's migration. 

1 Initial conditions 

The PPM code and initial conditions are very similar to those 
presented in Nelson et al. 1998. We begin with a one Mq pro- 
tostar fixed to the origin of our coordinate system. We assume 
that a disk of mass M^j = O.OSM© is contained between the 
inner and outer grid boundaries at 0.5 AU and 20 AU and that 
the disk is self gravitating. A second point mass (the 'planet') 
is set in a circular orbit at a radius 5.2 AU away from the pro- 
tostar and is free to migrate through the disk in response to 
gravitational forces. No other forces act on the planet and it 
does not accrete mass from the disk. In different simulations, 
we investigate migration rates of different planet masses. 

The disk mass is distributed on a 128x224 cylindrical (r, 0) 
grid with a surface density given by a power law, S(r) = 
Si {lAU /r)^ , where Ei is determined from the assumed disk 
mass and p = 3/2. We assume an initial temperature profile 
with a similar power law, T{r) = Ti (lACZ/r)', where the 
temperature at 1 AU is Ti = 250 K and q = 1/2. These initial 
conditions produce a radial profile for which the minimum 
Toomre Q (of ~ 5) is found near the outer disk edge. The 
profile exhibits a steep increase in the inner regions due to 
the increased effects of pressure on the orbital characteristics 
there. A single component isothermal gas equation of state is 
used to derive pressure at each point in the disk. 

Velocities are determined assuming initial rotational equi- 
librium. Radial velocities throughout are assumed equal to 
zero, while angular velocities are determined by balancing the 
gravitational, pressure and centrifugal forces in the disk. 

2 Migration rates with varying planet mass 

Using the initial conditions outlined above we have completed 
a series of simulations, varying the planet mass assumed for 
each simulation. We show the effect of a IMj mass planet on 
the disk in figure |^. Within a few hundred years, the planet 
raises very large amplitude spiral structures which lead (trail) 
the planet radially inside (outside) its orbit radius and cause a 
gap to form around the planet. These structures are similar to 
those in previous work (Bryden et al. 1999, Kley 1999) where 
the planet's trajectory was fixed to a single orbit radius. 

The orbital trajectory of the planet (figure ^) is strongly af- 
fected by the gravitational torques from the spiral structures. 
We show the After a ~ 50 — 100 yr period in which the planet 
first builds spiral structures, the migration rate is constant for 
the next 500 yr, then drops to zero and the planet reaches a 'fi- 
nal' orbit radius of about 3.5 AU, after 800 yr. The migration 
rates fit for simulations with different disk masses are shown 
in figure H, and are valid for theoeriod for which the migra- 
tion rate is constant (see figure Q), i.e. the period for which 
the migration can be considered 'Type F. The fitted migration 
rates are very rapid: about 1 AU per thousand years, so that 
migration in to the stellar surface would occur before the end 
of the disk lifetime of 10® yr. Still higher mass planets move 
so quickly through the disk that they 'outrun' their own gap 
formation efforts and fall into the star. With more realistic 
initial conditions (a 2Mj planet should already have a gap), 
we expect this phenomenon to go away. 




Figure 1 : The surface density of the disk after 775 yr of evolu- 
tion. As shown by the deep blue region in the plot, the planet 
has built a deep gap in the disk, but some spiral structure 
remains present in the gap. The planet has only just stopped 
its rapid migration through the disk at this time. 
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Figure 2: The orbital radius of the IMj planet shown in 
figure as a function of time. The planet begins at 5.2 AU 
and migrates inward to ~ 3.5 AU, before its inward trajectory 
is affected by the formation of a sufficiently deep and wide 
gap. 



Over the course of the first several hundred years of evolu- 
tion, planets with mass less than 2Mj hollow out a deep gap in 
the disk, which extends all the way around the star. In figure 
^ we show the azimuth averaged surface density structure at 
several points during the evolution of the system shown above. 
The gap forms quickly after the beginning of the simulation 
and within 500 years has hollowed out a region about 3 AU 
wide. The surface density near the planet (200-300 gm/cm^) 
is a factor ten below the initial profile and less (100 gm/cm^) 
at its deepest, just outside the planet's orbit. By the end of the 
simulation, the gap has deepened to a factor of 100 less than 
its unperturbed profile and continued to get deeper even at 
the end of our simulations. It does not substantially increase 
its width after initial formation however. 

The existence of a gap eventually causes the migration to 
slow and, if it becomes deep and wide enough, decrease to a 
rate defined by the viscosity of the disk ('Type II' migration) 
rather than to dynamical processes like gravitational torques 
('Type I' migration). From figure W we can determine the ap- 
proximate disk conditions which define the transition between 
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Figure 3; Migration rates of planets of varying mass. Points 
marked ' X ' are able to evacuate a gap large enough and deep 
enough to halt their migration. Planets greater than 1.5Mj 
'outrun' their own gap formation efforts and rapidly migrate 
inward to the inner grid boundary. 



Type I and Type II migration. From figure |2| we see that the 
migration rate decreases starting after about 5-600 yr of evo- 
lution. Comparison with figure W shows that the onset of the 
transition occurs when the gap is 3 AU wide and has surface 
density Sgap ~ 200 — 300gm/cm^. The transition is concluded 
and further rapid orbital decay is suppressed by the time the 
system has evolved for 700 yr, when the surface density is 
Sgap ~ 100 gm/cm^. 




Radius (AU) 



Figure 4: Azimuth averaged density structure at several times 
during the evolution. 

The values of surface density and gap width required for the 
onset and completion of the transition to Type II migration 
are typical of of each of the simulations we have performed, 
withouth regard to planet mass. Simulations with 0.3-0.5Mj 
planets are able to enter the transition to Type II migration 
but over the duration of our simulations (1800 yr) do not com- 
plete the transition. We continue to evolve these simulations 
further in time to determine their ultimate fate. 

3 The relative importance of the disk close to the 
planet 

Linear theory (e.g. Takeuchi et al. 1996) predicts that the 
most important Fourier components of the spiral patterns 
raised in the disk will be those with azimuthal wavenumber 
m=20— 40, corresponding to an azimuthal wavelength near the 
planet of about 1 AU. The Lindblad resonances of these pat- 
terns will be at a distance radially inward and outward from 
the planet of Rlr = a{l ± l/m)^/^ ^ 0.3 AU, where a is the 
semi-major axis of the planet. Another relevant parameter is 
the Hill radius, Rh = a{Mj/3M,y^^ ^ 0.3 AU, defining the 
sphere of influence of the planet. Further, the z structure of 
the disk becomes important on the same spatial scales because 
in 2D the gravity will be effectively 'amplified' by the assump- 
tion that all the disk matter is in the z = plane rather than 
some at high altitudes more distant from the planet. Unfor- 



tunately each of these values are similar to the grid resolution 
size scale of computationally affordable simulations. 

In order to characterize both physical and numerical effects 
we have performed a series of simulations in which we vary 
the effective gravitational softening length of the planet, ef- 
fectively 'turning on' or off the effect of matter very close to 
the planet. The migration rates obtained from such a series of 
simulations are shown in figure |5| We again assume the initial 
conditions as above, with a planet of mass 0.3Mj . 
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Figure 5: The migration rate of a 0.3Afj planet with vary- 
ing assumed gravitational softening length. The planet's Hill 
radius is shown with a vertical dashed line. From the plot 
it is clear that the most influenetial portion of the disk for 
migration is within 1-2 Hill radii of the planet. 

The migration rate increases by a factor of about five as the 
softening decreases from 0.5 AU to 0.1 AU, clearly showing 
the importance of the distribution of disk matter close to the 
planet. The largest increases occur as e decreases to the size 
of the Hill radius or smaller. Below e = 0.08 AU, or half the 
size of one grid cell, the migration rate is numerically unstable 
and slight changes in the softening change the migration rate 
from zero to 1 AU/lOO yr (below the bottom of the plot). The 
dependence of the migration on the disk matter very near the 
planet is in qualitative agreement with with the conclusion of 
Ward (1997), who showed that the most important contribu- 
tion to the migration will come from the disk matter within 
about one disk scale height (in our case ~ 0.3 AU near the 
planet) of the planet's radial position. These results show that 
very high resolution three dimensional simulations of the re- 
gion around the planet will be required in order to understand 
the migration rates of a planet through the disk. 

4 Summary 

In the course of this study we have 

• Shown that planets evolved in a 2D disk without a gap 
migrate inward through a substantial fraction of their initial 
semi-major axis radius on timescales ~ 1000 yr. 

• Shown that planets with masses higher than ^ 0.5Mj 
can open a gap sufficiently wide and deep to drastically slow 
their migration through the disk (i.e. transition to Type II 
migration). Conditions for beginning the transition are that 
surface density of the disk near the planet be about 200-200 
gm/cm^ and that the gap be ~ 3 AU wide. Conditions for 
completing the transition to Type II migration are surface 
densities near the disk of ~ 100 gm/cm^. 

• Demostrated the critical importance of very high spatial 
resolution of the disk near the planet required for correct evo- 
lution of the planet's migration. 
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